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Introduction 
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Background 


•  Transverse  jets  (jets-in-crossflow)  have  been 
studied  extensively  in  academia  and  industry 

•  Applications: 

•  Smokestack  dispersion 

•  Gas  turbine  burners 

•  VTOL  Aircraft 

•  Use  is  well  established  in  mixing  of  two 
dissimilar  fluids  for  combustion  devices 

•  Certain  regimes  of  multiple  jets  have  not 
been  studied  and  are  currently  of  interest  to 
Air  Force  research  goals 
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Problem  Statement 


Temperature  (°C) 

-g.GQae+Ui 
*8*150c+n3 
+0.500e+02 
-3+250e+02 
+  0 . 000e+02 

-7*5GOe+02 
-J*2Mlfi+fl2 
-  7 .  C00e  +  02 
-S*750e+C2 
+6*  SDOe+GZ 
--€.2026+02 


•  Temperature 


Accumulated  inelastic  strain 


6.961*- 03 
2.000e-02 
l,750e-03 
1 . 500e-03 

I. 2506-03 

i.o  x  *-■: 

7, £DOe-0fl 
5.0O0e-0fl 
2. 5006-04 

J. uixe-UL 


•  AF  seeking  to  develop  reusable  high-performance  liquid  rocket  engine 

•  Hindrances  to  reusability  are  turbine  cycle  fatigue  and  oxygen 
compatibility 

•  Turbine  cycle  fatigue  issues 

•  Density 


•  Species 
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Project  Objective 


•  Investigate  mixing  behavior  of  multiple  confined  JICF  to 
determine  local  optima 

•  Control  parameters  include  velocity  ratio,  diameter  ratio 
and  number  of  jets 

•  Non-dimensional  (ND)  parameters  of  interest  for  scaling 
include  momentum  flux  ratio,  momentum  ratio  and 
others 
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Confined  Jet  Configuration 


•  Studies  have  historically  focused  on 
gas  turbine  applications  involving 
rectangular  ducts  or  free  jets 

•  Holdeman  et  al.  studied  8-12 
cylindrical  orifaces  in  a  confined  can 

•  Downstream  profiles  similar  when 
scaled  with  orifice  spacing  (S/Dc)  and 
the  square  root  of  the  momentum  flux 
ratio 

5 

C  =  (— 

•  For  most  cases,  C  ~  2.5  (empirical) 

•  The  optimum  momentum  flux  ratio 
can  then  be  given  by: 

- ;(?)' 
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Unmixedness 


Gas  turbine  combustors  commonly  use 
temperature  as  a  means  of  analogizing 
concentration  since  they  are  both  scalar 
quantities.  The  non-dimensional  temperature 
can  be  defined  as: 


Using  optical  techniques,  concentration  can 
instead  be  analogized  by  light  intensity  (non- 
intrusive) 

Spatial  unmixedness  can  then  be  defined  from 
the  variance  and  normalized: 


U  = 


c2 


Where  m  is  the  pixel  count  and  C  is 
concentration  determined  from  calibration  data 
and  field  corrections 
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Relevant  Ratios 


Density  Ratio 

pj 

ni  = 

Pc 

Not  a  priority  in  water 
flow  experiment 

Velocity  Ratio 

•  Velocity  of  a  single  jet 
divided  by  crossflow 

Vj 

n2=  — 

vc 

Diameter  Ratio 

_  dJ 

•  Diameter  of  a  single  jet 

77-3  d 

divided  by  crossflow 

c 

Number  of  Jets 

TTq.  —  Tlj 
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Methodology 


•  Fix  jet  mass  flow  rates  by  cavitating  venturis 

•  Adjust  velocity  ratio  by  changing  core  mass  flow 
rate/core  velocity  for  given  geometry 


•  Fix  core  diameter  via  standard  test  section 

•  Adjust  jet  diameters  via  interchangable  inserts 


•  Maximum  number  of  jets  set  by  test  section 

•  Use  "blank"  inserts  to  reduce  the  number  of  jets 
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Design  Space 


•  How  to  determine  relative  ranges  of 
velocity  ratio,  diameter  ratio,  and  number 
of  jets? 

•  Holdeman  et  al  studied  numbers  of  jets 
from  8-16 

•  Knowledge  at  lower  n-values  important 
for  AF  research  goals 

•2,3,6  decided  upon  for  symmetric 
injection 

•  Previous  JICF  work  estimates  diameter 
ratios  between  0.1  -0.3 

•Used  knowledge  of  pre-existing  water  flow 
facility  to  approximately  size  core  diameter 
and  jet  diameters 


Jets  start  to 
impede  at 
D/Dc  =  0.833 
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Design  Space 


•  Previous  JICF  work  studies  J  values  from  5-200 

•  AF  research  focused  on  J  *  1 0-50 

•  Relevant  velocity  ratios  set  as  square  root  of 
these  values 


•Previous  JICF  work  estimates  diameter  ratios 
between  0.1 -0.3 

•  Used  knowledge  of  pre-existing  water  flow 
facility  to  approximately  size  core  diameter  and 
jet  diameters 

•  Based  on  mass  flow  available,  selected  0.12- 
0.21  as  possible  diameter  ratios 
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Diameter  Ratio 


Facility  Flow  Rate  Limit  (~2.3  GPM) 
Range  of  Investigation 


Njets  =  2 
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Design  Space 


^2  =  [VTo  V20  V30  V40 

7T3  =  [0.12  0.165  0.21] 

7T4  =  [2  3  6] 

•  Changing  diameter  ratio  or  number  of  jets  requires 
hardware  change 

•  Changing  number  of  jets  is  slightly  easier,  since  certain 
inserts  will  be  left  in  place  and  others  will  be  blocked  off  by 
using  “blanks” 

•  Diameter  ratio/number  of  jets  assigned  fewer  values  to 
minimize  hardware  changeover  while  still  covering 
parameters  of  interest 

•  Velocity  ratio  easiest  to  change,  assigned  more  test  points 
in  order  to  increase  strength  of  correlations 


v50  ] 

5  X  3  X  3  =  45  tests 


Distribution  A.  Approved  for  public  release;  distribution  unlimited 


15 


Apparatus 
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Water  Flow  Loop 


•  Pre-existing  hardware,  being  remodeled  to  accommodate  JICF 

•  Recirculating  loop,  tank  open  to  atm.,  manually  operated 

•  Approx,  capabilities:  100  psi  discharge,  100  gpm,  inlet  temp  10  -  40  °C 

•  Good  control  of  fluid  temp,  pressure  &  flow  rate 


Heat 

Exchanger 
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Total  Head  [psi] 


Water  Flow  Loop 


•  Water  Tank 

•  60  gal  cap 

•  Open  to  atm 

•  Pump 

•  Ebara  A3U32-200 

•  10  hp  (variable  RPM) 

•  Q@  10-120  GPM 

•  H  @  10-100  psi 


Operating  Curves  for  Ebara  Pump  Model  A3U32-200-10hp 
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Water  Flow  Loop 


•  Line  sizes  (D  =  1"  or  2")  allows  low  or  high  Q. 
operation  while  controlling  velocity 

•  Using  1"  line  for  most  tests 
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Modified  Flow  Loop 


•  Introduced  injection  piping 

•  Water  used  for  shakedown/start-up 


Heat 

Exchanger 


Fluorescein  used  for  data  collection 


Tank 


Waste 
Collection 


U 
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Turbulent  Run-Up 


•  Orifice  plate  between  inlet  and  turbulent  run  up  to 
increase  hydraulic  resistance  of  the  test  section 

•  1.5"  PVC  nominal 
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LDV  Window 


Provides  optical  access 
before  test  section 

LDV  crosshairs  will 
traverse  the  window 
laterally  to  measure  the 
inlet  velocity  profile 
(mean  and  turbulent 
fluctuations) 

Characterize  for  all 
potential  core  inflow 
velocities 
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Injection  Block 


•  Focal  point  of  jets 

•  Nylon,  machined  on-site  by 
AFRL techs 

•  Injection  cross  section  is 
same  as  turbulent  run-up 
(1.5"  nominal) 
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Jet  Inserts 


All  jet  inserts  machined  the  same  for  interchangability 
Characterize  using  LDV  system 

Ajayi,  Papadopoulos  and  Durst  (1998)  concluded  that  x/D  *  20  sufficient 
for  velocity  normalization  (5300  <  Re  <  10000) 


D/Dc 

L/Dj 

Rei 

0.12 

35.6 

10430 

0.165 

25.9 

7590 

0.21 

20.3 

5960 
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Injection  Manifold 


Pressure  drop  across  injection  manifold  is  minimal 
Variation  between  lines  is  <  0.07  psi 

Resulting  AQ  between  lines  is  <  1% 

■ 


Ball  valves 


■  ■ 


tsm  ■ 


To  injectors 


Cavitating 

venturis 


Tees  for 
pressure 
taps 


From 

injection 

tanks 
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Measurement  Window 


Laser 

Access 


Collection 

Plenum 


* 


Measurement 

Plane 


Plexi  Test 
Section 

nwm 

ML. 
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Data 

Camera 


Distribution  A.  Approved  for  public  release;  distribution  unlimited 


26 


Flow  Measurement 
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miniLDV  System 


Lenses,  beam  splitter,  and 
photodetector  housed  in  a 
single  unit 

Powered  by  60  mW  diode  laser 
(output  power  44  mW) 

Measurement  specs: 


•  1  mm/s  -  300  m/s 

•  Accuracy  99.7% 

•  Delivers  U  and  u' 


Image  Credit:  Measurement  Sciences 
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LDV  Profiles 
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LDV  Profiles 


Q-2.34GPM  ♦LDV  (Re  =  4468) 
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A  den  Toonder  LDV  (Re  =  4900) 

■ 

n  ft  ▲ 

U.O 

▲ 

n  c 

U.O 

n  a 

VJA 

n  9  A 

°-2  % 

o  , K 

0. 

U  i  i  i  I  w 

5  0.4  0.3  0.2  0.1  0 

-0.2 

Distribution  A.  Approved  for  public  release;  distribution  unlimited 


30 


41* 


LDV  Profiles 
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Planar  Laser  Induced  Fluorescence 


Laser  spread  to  a  sheet  using 
a  cylindrical  lens 

Laser  energy  excites 
fluorescent  agent  seeded 
into  fluid  of  interest 

Calibrate  using: 

•  Laser  Sheet  Intensity 

•  Injection  fluid 
concentration 

•  Exposure  time 


Image  Credit:  R.  McLaughlin,  University  of  North  Carolina,  Chapel  Hill 
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Planar  Laser  Induced  Fluorescence 


•  Fluorescent  agent  is  sodium 
fluorescein  (C2oH10Na205) 

•  Benefits: 

•  High  quantum  yield  (~94%) 

•  Water  soluable 

•  Relatively  safe 

•  Other  chemicals  considered  were 
Rhodamine  (G  and  6B); 
downselected  for  safety  and  ease 
of  handling 

•  Absorption  spectrum  overlaps 
sufficiently  with  Ar-lon  laser  @  514 
nm 


Absorption  Max 

~494  nm 

Emission  Max 

~521  nm 

Yellow  filter  eliminates  these 
bands  from  data  collection 


-r  i  i  i  *  ii - 1 - -4- 

350  400  450  500  550  600  650 

Wavelength  (nm) 


Image  Credit:  Invitrogen 
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Fluorescence  emission 


In-Situ  Calibration 


Raw  field  intensity 

•  Data  set  taken  per  test 

•  3-4  images  per  test 
Dark  field  intensity 

•  Taken  before  each  test 

•  Background  correction  for 
camera  pixels 

Reference  field  intensity 

•  Taken  daily 

•  Fluid-filled  test  section  at  pure 
jet  fluid  concentration 

•  Accounts  for  attenuation,  beam 
shape  and  other  environmental 
factors 


j  _  I  RAW  IDF 

*NQRM  ~  ~T  ~T 

*REF  ~  lDF 


Sample  reference  image 
Highpass  Filtered  at  530  nm 
C  =  4.8xl0'8  mol/L 
f-stop  =  11 
Exp  time  =  2  sec 
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Results 

We’re  not  kidding! 
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Proposed  Results 


Holdeman  deviation 


—  rij  <  n 


■critical 


(=6) 


1  fCntf 
fopt  2  \  n  / 


A  change  in  the  optimum 
unmixedness  with  the  ND 
parameter  of  interest  would 
indicate  poor  scaling 


Momentum  Flux  vs  Number  of  Jets 

EXAMPLE  ONLY. 


Optimum 
Momentum 
Flux  Ratio,  Jopt 


—♦—Holdeman 

-B-Dj/Dc  =  0.12 

/ 

Dj/Dc  =  0.165 

-•-Dj/Dc  =  0.21 

Show  how  unmixedness 
varies  with  ND  parameters 


Concentration  images 

Analysis  of  macro-level  spatial 
flow  features 


20  40 

ND  Parameter  of  Interest 
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Un  mixed  ness  Urimixedness 


Unmixedness 
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Unmixedness 


Unmixedness 
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Unmixedness 


Unmixedness 


Mixing  vs  J 
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Sample  Case  -  N1_J10_D12 


Processed  image  with  contours 


Raw  image 
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Sample  Case  -  N6_J30_D12 
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Holdeman  Plot 


Mixing  vs  J 
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Conclusions 


*  Sparse  multi-JICF  behave  much  how  we  expect 

*  Break-off  point  with  Holdeman  scaling  slightly 
higher  than  theory  predicts  (2tt  ~  6.28) 

*  Diminishing  returns  moving  from  3  to  6  jets 

*  Take-away:  there  is  an  optimum  J-value  for  the  six- 
jet  configuration  near  J  =  30 
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Future  Work 


*  Obtain  remainder  of  data  for  higher  diameter  ratios* 

*  Look  at  other  axial  locations  for  promising  cases* 

*  Re-evaluate  data  anomalies 

*  Write  a  thesis?  (no  big  deal) 


*Themis  work 
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Lessons  Learned  (Gripes) 


*  Government  work  is  just  as  slow  as  everyone  jokes 

*  Experimentation  is  full  of  headaches 

—  Tight  fitting  of  jet  inserts  with  injection  block 
—  Determining  appropriate  concentration 
—  Rusty  tanks 
—  Old  pressure  transducers 
—  Leaks,  leaks,  leaks 

—  Having  to  explain  to  safety  that  fluorescein  is  less  dangerous  than 
vacuum  pockets 

—  Fixing  Labview  (and  logistics  of  dealing  with  IT) 

—  Doing  it  all  on  a  shoestring  budget 

*  Just  use  gases,  they're  way  more  simple 
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Project  Themis: 
Water  Visualization  Study 


Cal  Poly 

Dr.  Tina  Jameson  (chair) 
Dr.  Dianne  DeTurris 
Dr.  David  Marshall 


AFRL 

Dr.  Rich  Cohn 
Nils  Sedano 
Dr.  David  Forliti 
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Back-Ups 
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Transverse  Jets 


Unconfined  single  transverse  jets  (single  phase) 

•  Exhaust  stacks,  V/STOL  aircraft,  3D  canonical  flow 

•  Studied  extensively  in  literature 


From  Fric  and  Roshko 


Jet  ihcar'byrr 
vorlkxi 


•  Fric  and  Roshko  (1 994)  cited  358  times 

•  Smith  and  Mungal  (1998)  cited  201  times 

•  Yuan,  Street  and  Ferziger  (1 999)  cited  1 55  times 

•  Margason  review  article  (1993)  cited  165  times 

Multiple  vortex  structures 
Scaling  laws 


Unconfined  single  transverse  jets  (two  phase) 

•  Atomization  of  liquid  jets 

•  Supersonic  gas-phase  velocities  (Scramjet  application) 


Confined  multiple  transverse  jets 

•  Motivated  to  support  gas  turbine  combustor  design 

•  Holdeman  (NASA  Glenn)  and  coworkers  major 
contributor 
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Transverse  Jet  Anatomy 


Counter-rotating  Vortex  Pair 
(CVP) 

Horseshoe  vortex  forms  at 
forward  stagnation  point 

Wake  vortices  shed  from  jet 

Occur  over  a  wide  range  of 
Reynolds  numbers 


Counter- rota  ri  rig 

vortex  pair 


Wake  vuiliuci 


Wal  l 


Jet  shciir-liiyiT 
vortices 


Image  Credit:  Fric  and  Roshko  (1994) 
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Problem  Statement 


Current  US  Launch  Systems 


Delta  IV  (Boeing/ULA) 

•  First  Stage:  1  RS-68 

•  Propellants:  LOX/LH2 


Atlas  V  (Lockheed  Martin/ULA) 

•  First  Stage:  1  RD-180 

•  Propellants:  LOX/RP-1 
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Parameters  and  Scaling  Laws 


Physical  argument  for  the  Holdeman  parameter 
For  good  mixing,  jet  penetration  ~  Radius 


l2U  2  =  D  2U  2 

Ctc  c  j  j 


42  %Dj ) 


blockage 


Djr  ccRc 


Jet  length  scale  is  rDj: 

•Jet  strongly  deflected  in  a  streamwise  length  of  rDj 
•Jet  growth  to  a  width  of  ~rDj 

Good  mixing: 

Jets  far  enough  apart  to  allow  CVP  to  form 

Jets  close  enough  to  begin  to  merge  when  they  approach  a  size  of  rDj 


DjV  x  S 
*  Djr2  oc  Rc 


Sr  oc  Rc  argue  r  oc  J 


1/2 
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ND  Groups 


•  Used  Buckingham-Pi  to  determine  relevant 
dimensionless  parameters 

•  Relevant  equations: 

Continuity  Momentum  Total  Momentum 

rhj  =  PjAjVj  pjAjVj2  PjAjVj2rij 

Pc^-c^c  Pc^-c^c  Pc^-c^c 

•  Recurring  Variables: 

Pj  Pc  Vj  Tty 

•  Other  Variables 


Pj  Pc 
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Parameters  of  Interest 


ND  Parameter 

Equations 

Notes 

Momentum  Flux  Ratio 

Prominent  in  literature 

"Kappa" 

K  = 

ll\2  ( l') 

vj  \dc) 

Proportional  to  rD  scaling 

"Lambda" 

X  =  | 

my 

Proportional  to  rD  scaling 

"Alpha" 

a  = 

m 

\pc/\] 

vj  \dj 

Proportional  to  r2D  scaling 

"Epsilon" 

e  =  ( 

M  (vj_ 

Pc' 

my 

Proportional  to  r2D  scaling 

Momentum  Ratio 

H 

fp£\  p 

<pc)  \vc 

mi 

Prominent  for  certain  jet  types 

Total  Momentum  Ratio 

i />  =  ( ‘ 

0c)\vc> 

my 

Accounts  for  jet  interaction 
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Laser  Doppler  Velocimetry  (LDV) 


Use  interference  patterns  and 
Mie  scattering  to  calculate  bulk 
flow  velocity 

Using  geometry,  the  fringe 
spacing  6  can  be  determined 
from  laser  wavelength  and 
convergence  angle  0: 

X 

5  =  — w 


2  sin  (2) 


Whenever  a  scattering 
particle  crosses  one  of  these 
interference  planes,  it  will  be 
sensed  by  a  photodetector 


Image  Credit:  eFunda 
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Laser  Doppler  Velocimetry 


•  Particles  will  periodically  pass  through  interference  planes  (not  continuously) 

vN 

•  This  corresponds  to  the  Doppler  burst  frequency:  ^D=~s  where  vN  represents 
the  component  of  velocity  normal  to  the  interference  plane 


•  Limitations: 

•  Cannot  predict  direction 

•  True  velocity  must  be  estimated 
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LDV  Profiles 


n  9c 

0.2 

0.15 

0.1 

0.05 

ft 

u=  l.bl  or  M 

-r 

♦  LDV  (Re  =  3000) 

j  Evan  Doorne 

u 

1  I  T  i 

<►  * 

■  1 

M 

1  1  T  1.* 

.ill 

B1 

u 

c 

i  i  i  i  ▼  i 

)  0.1  0.2  0.3  0.4  0.5  0.6 
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LDV  Profiles 


0  0.1  0.2  0.3  0.4  0.5  0.6 
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LDV  Profiles 
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Sample  PLIF  Setup 


*  / 


Bl- 


Laser 

(Ar-lon) 


Test  Article 


Redirection 

Mirrors 


t  i  ^ 


i.r 


CCD 

Camera 


I#  ^ 


Yellow  Filter 
(-530  nm) 
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General  PLIF  Considerations 


.  Mu,dpH>9 

•  Photobleaching 

•  Fluid  purity 

•  Laser  power 


beam/sneet  attenuation 

Camera 

Environment  and 
background  correction 


Laser  Attenuation  vs.  Distance 


♦  100  mW 
■  1 50  mW 

200  mW 
250  mW 
x  300  mW 

•  350  mW 
+  400  mW 
=  450  mW 


Shadow  streaking  as  a  result  of  fluid  impurities 


100  200 

Distance  (pixels) 


400 
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